Quantum Error Correction of Time- Correlated Errors 



Feng Lu and Dan C. Marinescu 
School of Electrical Engineering and Computer Science 
. . . _ University of Central Florida, P. O. Box 162362 

^ ■ Orlando, FL 32816-2362 

. Email: (Ivfeng, dcm)@cs. ucf.edu 

<N : 

^ . February 9, 2008 

' Abstract 

The complexity of the error correction circuitry forces us to design quantum error correction 
' codes capable of correcting a single error per error correction cycle. Yet, time-correlated error 

^ I are common for physical implementations of quantum systems; an error corrected during the 

■ previous cycle may reoccur later due to physical processes specific for each physical implemen- 
(N ■ tation of the qubits. In this paper we study quantum error correction for a restricted class of 

time-correlated errors in a spin-boson model. The algorithm we propose allows the correction 

■ of two errors per error correction cycle, provided that one of them is time-correlated. The al- 
' gorithm can be applied to any stabilizer code when the two logical qubits | Ol) and | 1l) are 
, entangled states of 2" basis states in H2" . 

^ : 

1 Quantum Error Correction 

■ 

. Quantum states are subject to decoherence and the question whether a rehable quantum computer 

^ ' could be built was asked early on. A "pure state" | (/s) = ao | 0) -|- ai | 1) may be transformed as a 

result of the interaction with the environment into a "mixed state" with density matrix: 

p=|aonO)(0| + |aini)(l|. 
^ ■ Other forms of decoherence, e.g. leakage may affect the state probability amplitude as well. 



The initial thought was that a quantum computation could only be carried out successfully if its 
duration is shorter than the decoherence time of the quantum computer. As we shall see in Section 
m the decoherence time ranges from about 10^ seconds for the nuclear spin embodiment of a qubit, 
to 10~^ seconds for quantum dots based upon charge. Thus, it seemed very problematic that a 
quantum computer could be built unless we have a mechanism to deal periodically with errors. Now 
we know |16j that quantum error correcting codes can be used to ensure fault-tolerant quantum 
computing; quantum error correction allows us to deal algorithmically with decoherence. There is 
a significant price to pay to achieve fault-tolerance through error correction: the number of qubits 
required to correct errors could be several orders of magnitude larger than the number of "useful" 
qubits 0. 

We describe the effect of the environment upon a qubit as a transformation given by Pauli 
operators: (i) the state of the qubit is unchanged if we apply the aj operator; (ii) a bit-flip error is 
the result of applying the transformation given by ; (hi) a phase-flip error is the result of applying 
the transformation given by cr^; and (iv) a bit-and-phase flip error is the result of applying the 
transformation given by ay = iUxOz- 

A quantum error correcting scheme takes advantage of entanglement in two ways: 

• We entangle one qubit carrying information with {n— 1) other qubits initially in state | 0) and 
create an n-qubit quantum codeword which is more resilient to errors. 
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• We entangle the n qubits of the quantum codeword with ancilla qubits in such a way that 
we can measure the ancilla qubits to determine the error syndrome without altering the state 
of the n-qubit codeword, by performing a so called non demolition measurement. The error 
syndrome tells if the individual qubits of the codeword have been affected by errors as well as 
the type of error. 

• Finally, we correct the error(s). 

Even though the no-cloning theorem prohibits the replication of a quantum state, we are able 
to encode a single logical qubit as multiple physical qubits and thus we can correct quantum errors 
[Is] . For example, we can encode the state of a qubit 

\^) = an\0)+ai I 1} 

as a linear combination of | 00000} and | 11111): 

\ip) ^ao\ 00000) + ai | 11111). 

Alternatively, we can encode the qubit in state | ip) as: 

\'p)=ao\ Oi) +ai I U), 

with I Ol) and | 1^) expressed as a superposition of codewords of a classical linear code. In this 
case all codewords are superpositions of vectors in a Hilbert space of dimension 2^. Steane's 
seven-qubit code |17j and Shor's nine qubit code [3] are based upon this scheme. 

When we use the first encoding scheme, a random error can cause departures from the subspace 
spanned by | 00000) and | 11111). We should be able to correct small bit-flip errors because the 
component which was | 00000) is likely to remain in a sub-space Co C spanned by the six vectors: 

I 00000), I 00001), I 00010), I 00100), I 01000), I 10000) 

while the component which was | 11111) is likely to remain in a sub-space Ci C 7^2= spanned by the 
six vectors, 

I 11111), I 11110), I 11101), I 11011), I 10111), I 01111). 
The two subspaces are disjoint: 

Co n Ci = 

thus we are able to correct a bit- flip error of any single physical qubit. This procedure is reminiscent 
of the basic idea of classical error correction when we determine the Hamming sphere an n-tuple 
belongs to, and then correct it as the codeword at the center of the sphere. 

A quantum code encodes one logical qubit into n physical qubits in the Hilbert space 7^2" with 
basis vectors {| 0), | 1), . . . | i), . . . | 2" — 1)}. The two logical qubits are entangled states in 7^2": 

|Ot)=^a, |i) and |1l)=5IAI*) 

i i 

A quantum error correcting code must map coherently the two-dimensional Hilbert space spanned 
by I Ol) and | 1l) into two-dimensional Hilbert spaces corresponding to bit-flip, phase-flip, as well 
as bit-and-phase flip of each of the n qubits to ensure that the code is capable of correcting the three 
types of error for each qubit. The quantum Hamming bound: 

2{3n + l) < 2". 

established by Laflamme, Miquel, Paz, and Zurek in [TT] allows us to say that n = 5 is the smallest 
number of qubits required to encode the two superposition states | Ol) and | 1^), and then be able 
to recover them regardless of the qubit in error and the type of errors. 
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The same paper introduces a family of 5-qubit quantum error correcting codes. The code Q 
with the logical codewords 

Ol) = i (I 00000)+ I 10010}+ I 01001}+ I 10100}+ I 01010)- I 11011)- I 00110}- I 11000} 

- I 11101}- I 00011}- I 11110)- I 01111)- I 10001)- I 01100}- I 10111}+ I 00101)) 

(I 11111)+ I 01101}+ I 10110}+ I 01011}+ I 10101)- I 00100)- I 11001}- I 00111} 

- 1 00010}- 1 11100}- 1 00001)- 1 10000)- I OHIO)- | looii}- | oiooo}+ 1 iioio)) 

is a member of this family. The code Q is a perfect quantum error correcting code because a logical 
codeword consists of the smallest number n of qubits to satisfy the inequality from [11| . Recall that 
a perfect linear code [n, fc, d] with d = 2e + 1 is one where Hamming spheres of radius e are all disjoint 
and exhaust the entire space of n-tuples. 

Different physical implementations reveal that the interactions of the qubits with the environment 
are more complex and force us to consider spatially- as well as, time-correlated errors. If the qubits 
on an n qubit register are confined to a 3D structure, an error affecting one qubit will propagate to 
the qubits in a volume centered around the qubit in error. Spatially-correlated errors and means 
to deal with the spatial noise are analyzed in [J |31 [^. An error affecting qubit i of an n-qubit 
register may affect other qubits of the register. An error affecting qubit i at time t and corrected 
at time t + A may have further effect either on qubit i or on other qubits of the register. The error 
model considered in this paper is based upon a recent study which addresses the problem of reliable 
quantum computing using solid state systems |13j . 

There are two obvious approaches to deal with quantum correlated errors : 

• (i) design a code capable to correct these two or more errors, or 

• (ii) use the classical information regarding past errors and quantum error correcting codes 
capable to correct a single error. 

When a quantum error correcting code uses a large number of qubits and quantum gates it 
becomes increasingly more difficult to carry out the encoding and syndrome measurements during 
a single quantum error correction cycle. For example, if we encode a logical qubit using Shor's nine 
qubit code we need 9 physical qubits. If we use a two-level convolutional code based upon Shor's 
code then we need 81 physical qubits and 63 = 7 x 9 ancilla qubits to ensure that the circuit for 
syndrome calculation is fault-tolerant. While constructing quantum codes capable of correcting a 
larger number of errors is possible, we believe that the price to pay, the increase in circuit complexity 
makes this solution undesirable; this motivates our interest in the second approach. 

2 Time-Correlated Quantum Errors 

Figure [T] illustrates the evolution in time of two qubits, i and j for two error correction cycles. The 
first error correction cycle ends at time t2 and the second at time t^. At time ti qubit i is affected by 
decoherence and flips; at time ^2 it is flipped back to its original state during the first error correction 
step; at time ts qubit j is affected by decoherence and it is flipped; at time t^ the correlation effect 
discussed in this section affects qubit i and flips it to an error state. If an algorithm is capable to 
correct one correlated error in addition to a "new" error, then during the second error correction 
the errors affecting qubits i and j are corrected at time t^ . 

The quantum computer and the environment are entangled during the quantum computation. 
When we measure the state of the quantum computer this entanglement is translated into a proba- 
bility e that the measured state differs from the expected one. This probability of error determines 
the actual number of errors a quantum code is expected to correct. 

If Tgate is the time required for a single gate operation and Tdch is the decoherence time of a 
qubit, then Ugates-, the number of gates that can be traversed by a register before it is affected by 
decoherence is given by: 
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Figure 1: Time-correlated quantum errors. Two consecutive error correction cycles occurring at 
time t2 and ^5 are shown. At time a code designed to correct a single error will fail while a 
code capable to handle time-correlated errors will correct the "new" error of qubit j and the "old" 
correlated error of qubit i. 



Table 1: The time required for a single gate operation, Tgate, the decoherence time of a qubit, Tdch, 
and the number of gates that can be traversed before a register of qubits is affected by decoherence, 

'Agates • 



Qubit implementation 


'^dch{sec 


'^gate{sec) 


Agates 


Nuclear spin 


10* 


10- 


-3 


10^ 


Trapped Indium ion 


10-1 


10- 


-14 


1013 


Quantum dots/charge 


10-9 


10- 


-12 


103 


Quantum dots/spin 


10-6 


10- 


-9 


103 


Optical cavity 


10-5 


10- 


-14 


10^ 



_ Tdch 
'Agates — . 

Tgate 

Quantum error correction is intimately related to the physical processes which cause decoherence. 
Table [1] presents sample values of the time required for a single gate operation Tgate, the decoherence 
time of a qubit, Tdch, and Ugates, the number of gates that can be traversed before a register of qubits 
is affected by decoherence, for several qubit implementations [HI HH UHl . We notice a fair range of 
values for the number of quantum gate operations that can be performed before decoherence affects 
the state. 

The information in Table [U in particular the decoherence time, can be used to determine the 
length of an error correction cycle, the time elapsed between two consecutive error correction steps. 
The number of quantum gate operations limits the complexity of the quantum circuit required for 
quantum error correction. 

The Quantum Error Correction theory is based upon the assumption that the quantum system 
has a constant error rate e. This implies that once we correct an error at time tc, the system behavior 
at time t > is decoupled from events prior to t^. 

The Markovian error model is not consistent with some physical implementations of qubits; for 
example, in a recent paper Novais and Baranger [13] discuss the decoherence in a spin-boson model 
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which is apphcable, for instance, to quantum dots. The authors assume that the qubits are perfect, 
thus eventual errors are due to dephasing and consider a Unear couphng of the qubits to an ohmic 
bath. Using this model the authors analyze the threc-qubit Steane's code. They calculate the 
probability of having errors in quantum error correction cycles starting at times ti and t2 and show 
that the probability of errors consists of two terms; the first is the uncorrelated probability and the 
second is the contribution due to correlation between errors in different cycles (A is the period of 
the error correcting cycle): 



'2' 8(fl-<2)4 

We see that correlations in the quantum system decay algebraically in time, and the latest error will 
re-influence the system with a much higher probability than others. 

Only phase-flip errors are discussed in detail in [13] . The approach introduced by the authors is 
very general and can be extended to include other types of errors. Nevertheless, for other physical 
models, one does expect major departures from the results presented in |13j . 

3 Stabilizer Codes 

The stabilizer formalism is a succinct manner of describing a quantum error correcting code by a 
set of quantum operators [7]. We first review several concepts and properties of stabilizer codes. 

The 1-qubit Pauli group Gi consists of the Pauli operators, u/, ax-, Uy, and cr^ together with the 
multiplicative factors ±1 and ±i: 

Qi EE {±CT/, ±icr/, zbcTa;, iio-x, ±o-y, ±i(Ty, iCT^, ibitT^}. 
The generators of Qi are: 

The n-qubit Pauli group Qn consists of the 4" tensor products of cr/, Cy, and Gz with an overall 
phase of ±1 or ±i. Elements of the group can be used to describe the error operators applied to an 
n-qubit register. The weight of such an operator in Qn is equal to the number of tensor factors which 
are not equal to ct/. The stabilizer S of code Q is a subgroup of the n-qubit Pauli group, S C Qn- 
The generators of the subgroup S are: 

M = {Mi,M2...MJ. 

The eigenvectors of the generators {Mi, M2 . . . M^} have special properties: those corresponding to 
eigenvalues of -1-1 are the codewords of Q and those corresponding to eigenvalues of —1 are codewords 
affected by errors. If a vector \ ipi) £ W„ satisfies, 

M, I V^,) = (+1) I V^) VMj e M 

then I ipi) is a codeword, \ ipi) (z Q. This justifies the name given to the set S, any operator in S 
stabilizes a codeword, leaving the state of a codeword unchanged. On the other hand if: 

Mj I ipk) - (-1) I ipk). 

then I (pk) = | ipk), the state \ ipk) is a codeword \ ipk) Q affected by error E^. The error 
operators affecting codewords in Q, E = {Ei,E2...}, are also a subgroup of the n-qubit Pauli 
group 

ECQn- 

Each error operator E^ is a tensor product of n Pauli matrices. Its weight is equal to the number of 
errors affecting a quantum work, thus the number of Pauli matrices other than aj. 
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The coding space: 



g - {| e Hn such that Mj I = (+1) I ^,), VAfj e 5 } 
is the space of all vectors | ipi) fixed by S 

It is easy to prove that S* is a stabilizer of a non-trivial Hilbert subspace V2" C if E^nd only 

if: 

1.5'= {Si, S2, . . .} is an Abelian group: 

2. The identity matrix multiplied by —1 is not in S: 

-(af") ^ 5. 

If E is an error operator, and E anti-commutes with some element M G 5, then E can be 
detected, since for any | 1/7^) G Q: 

ME I V») = -EM I V,) = -E I 7/.,) 

Definition. The normalizer N{S) consists of operators E G 5n such that ES^E^^ G 5, VS^ G S. The 
distance d of a stabilizer code is the minimum weight of an element in N(S) — S with N{S) the 
normalizer of S. An [n, k, d] stabilizer code is an [n, k] code stabilized by S and with distance d; n 
is the length of a codeword, and k the number of information symbols. 

An [n, k, d] stabilizer code with the minimum distance d = 2e + 1 can correct at most e arbitrary 
quantum errors or 2e errors whose location is well known. In [5] this property of a stabilizer code is 
used for syndrome decoding for optical cluster state quantum computation. 

Given an [n, k, d] stabilizer code the cardinalities of the stabilizer S and of its generator M are: 



ls'l=2"-^ |M|=7i-fc 

The error syndrome corresponding to the stabilizer Mj is a function of the error operator, E, defined 
as: 

/m,(e) .g^Z2 /m,(e) - I ^ E} = 0, 

where [Mj , E] is the commutator and {Mj, E} the anti-commutator of operators Mj and E. Let 
/(E) be the (n — fc)-bit integer given by the binary vector: 



/(E) = (/mi(E)/m,(E).../m„_,(E)). 
This (n — A:)-bit integer is called the syndrome of error E. 

Proposition. The error syndrome uniquely identifies the qubit(s) in error if and only if the subsets 
of the stabilizer group which anti-commute with the error operators are distinct. 

An error can be identified and corrected only if it can be distinguished from any other error in 
the error set. Let Q{S) be the stabilizer code with stabilizer S. The Correctable Set of Errors for 
Q(S) includes all errors which can be detected by S and have distinct error syndromes. 
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Table 2: Single error operators for the 5-qubit code and the generator(s) that anti-commute with 
each operator. Subscripts indicate positions of errors, e.g. X3 means a bit-flip error on the third 
qubit, X3 = (7/ ® (7/ (g) fj^c (g) tj/ (X) CT/ . 



Error operator 


Generator(s) 


Error operator 


Generator(s) 


Error operator 


Generator(s) 


Xi 


M4 


Zi 


Ml, M3 


Yi 


Ml, M3, M4 


X2 


Ml 


Z2 


M2, M4 


Y2 


Ml, M2, M4 


X3 


Ml, M2 


Z3 


M3 


Y3 


Ml, M2, M3 


X4 


M2, M3 


Z4 


Ml, M4 


Y4 


Ml, M2, M3, M4 


X5 


M3, M4 


Z5 


M2 


Y5 


M2, M3, M4 



Corolary. Given a quantum error correcting code Q capable to correct e„ errors, the syndrome does 
not allow us to distinguish the case when more than e„ qubits are in error. If we can identify the Cc 
correlated errors in sysntem, the code is capable of correcting these e„ + Cc errors. 

Proof: Assume that -Fi, F2 cause at most e„ qubits to be in error, thus i^i, F2 are included in the 
Correctable Set of Errors of Q. Assuming errors -Fi and F2 are distinguishable, there must exist 
some operator M e S which commutes with one of them, and anti-conimutes with the other: 

Ff F2M = -MF[F2 

If we know the exact correlated errors E in the system, then: 

{E^Fi)'^{E^F2)M = {F[EE'^F2)M = F^ F2M = -MF^F2 = -~ M {E^ Fi)'^ {E^ F2) 
which means that the stabilizer AI commutes with one of the two errors E'^Fi, E'^F2 and anti- 
commutes with the other. So error E^ Fi is distinguishable from error E'^F2. Therefore, if we know 
the exact prior errors E, we can identify and correct any E^ Fi errors with the weight of Fi equal 
or less than 

For example, consider a 5-qubit quantum error-correcting code Q with n = 5 and fc = 1 from 
[TT] discussed in Section [TJ The stabilizer S of this code is described by a group of 4 generators: 

M = {Ml, M2, M3, M2} with the generators: 

Ml = CTa; (g) Cr^ (g) (7z (g) (Ta; (g) (T/ , M2 = CT/ (g) (Ta; (g) CTz (g) tJz ® (Ta, , 
M3 = (Ta; (g) (7/ ® (Ta; (g) Cr^ (g) , M4 = (T^ (g) (Ta; (g) (7/ (g) (Ta; ^ . 

It is easy to see that two codewords are eigenvectors of the stabilizers with an eigenvalue of (+1): 

I Ol) = (+1) I Ql) and M, | 1^) = (+1) | 1^), 1 < j < 4. 
Note also that M is an Abelian subgroup, each generator commutes with all the others. 

Table [2] lists single error operators and the generator(s) which anti-commute with each operator. 
For example, Xi anti-commutes with M4, thus a bit-flip on the first qubit can be detected; Zi 
anti-commutes with Mi and M3, thus a phase-flip of the first qubit can also be detected. Since each 
of these 15 errors anti-commute with distinct subsets of S we can distinguish individual errors and 
then correct them. An example shows that the code cannot detect two qubit errors; indeed, the two 
bit-flip error 

X1X2 = CTx ® CTx ® cri ® ai ® ai 

is indistinguishable from Z4 — (Ji®ai®ai®(Jz® o'l because both X1X2 and Z4 anti-commute 
with the same subset of stabilizers, {Mi,M4}, and give the same error syndrome. Therefore, the 
5-qubit code can correct any single qubit error, but cannot correct two qubit errors. 
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If we know the correlated error in the system, for example, X3, then, from Tabled it is easy to 
see that the errors XsE.i have distinct error syndromes for 1 < z < 15. Therefore we can identify 
these errors and correct them. 

4 Time-Correlated Errors and Stabilizer Codes 

Classical, as well as quantum error correction schemes allow us to construct codes with a well-defined 
error correction capability. If a code is designed to correct e errors, it will fail whenever more than 
e errors occur. 

From previous section we know that if time-correlated errors are present in the system, we can 
not always detect them through the calculation of the syndrome. The same syndrome could signal 
the presence of a single error, two, or more errors, as shown by our example in Section [3l we can 
remove this ambiguity when there are two errors and one of them is a time-correlated error. 

In this section we extend the error correction capabilities of any stabilizer code designed to 
correct a single error (bit-flip, phase-flip error, or bit-and-phase flip) and allow the code to correct 
an additional time-correlated error. The standard assumptions for quantum error correction are: 

• Quantum gates are perfect and operate much faster than the characteristic response of the 
environment; 

• The states of the computer can be prepared with no errors. 
We also assume that: 

• There is no spatial-correlation of errors, a qubit in error does not influence its neighbors; 

• In each error correcting cycle, in addition to a new error Ea that occurs with a constant 
probability Sa, a time-correlated error Eh may occur with probability eb{t). As correlations 
decay in time, the qubit affected by error during the previous error correction cycle, has the 
highest probability to relapse. 

Quantum error correction requires non- demolition measurements of the error syndrome in order 
to preserve the state of the physical qubits. In other words, a measurement of the probe (the ancilla 
qubits) should not influence the free motion of the signal system. The syndrome has to identify 
precisely the qubit(s) in error and the type of error(s). Thus, the qubits of the syndrome are either 
in the state with density | 0)(0 | or in the state with density | 1)(1 |, which represent classical 
information. 

A quantum non-demolition measurement allows us to construct the error syndrome Ticurrent- 
After examining the syndrome "Scurrent an error correcting algorithm should be able to decide 
whether: 

1. No error has occurred; no action should be taken; 

2. One "new" error, Ea-, has occurred; then we apply the corresponding Pauli transformation to 
the qubit in error; 

3. Two or more errors have occurred. There are two distinct possibilities: (a) We have a "new" 
error as well as an "old" one, the time-correlated error; (b) There are two or more "new" 
errors. A quantum error correcting code capable of correcting a single error will fail in both 
cases. 

It is rather hard to distinguish the last two possibilities. For perfect codes, the syndrome Sab 
corresponding to two errors, Ea and Ei,, is always identical to the syndrome Sc for some single 
error E^- Thus, for perfect quantum codes the stabilizer formalism does not allow us to distinguish 
two errors from a single one; for a non-perfect code it is sometimes possible to distinguish the two 
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syndromes and then using the knowledge regarding the time-correlated error it may be possible to 
correct both the "old" and the "new" error. 

We now describe an algorithm based on the stabilizer formalism capable to handle case 3(a) for 
perfect as well as non-perfect quantum codes. We assume that the quantum code has a codeword 
consisting of n qubits and uses k ancilla qubits for syndrome measurement. 

The outline of our algorithm: 

1. At the end of an error correction cycle entangle the qubit affected by the error with two 
additional ancilla qubits. Thus, if the qubit relapses, the time-correlated error will propagate 
to the two additional ancilla qubits. 

2. Carry out an extended syndrome measurement of the codeword together with the two additional 
ancilla qubits. This syndrome measurement should not alter the state of the codeword and 
keep the entanglement between the codeword and the two additional ancilla qubits intact. 

3. Disentangle the two additional ancilla qubits from the codeword and then measure the two 
additional ancilla qubits. 

4. Carry out the error correction according to the outcomes of Steps 2 and 3. 

Next we use the Steane code to illustrate the details of the algorithm. The generators of the 
Steane 7-qubit code are: 

M.I ^ aj ® ai ® ai ® ax ® (Tx ® Ux ® (Tx, M2 = ai ® ax ® (Jx ® (^i ® (^i ® ctx cTx, 
M3 — ax ® ai ® ax (ii ai ® ax ® ai ® ax, M4 — ai ® ai ® ai ® a^ ® az ® az ® az, 
M5 ~ai®az®az®ai®aj®az®az, Mg ~az®ai®az®aj®az®ai®az- 

1. Entangle the two additional ancilla qubits with the original code word. The time- 
correlated error may reoccur in different physical systems differently: a bit-flip could lead to a 
bit-flip, a phase-flip, or a bit-phase- flip. Therefore, we need to "backup" the qubit corrected during 
the previous error correction cycle in both X and Z basis. We entangle the qubit affected by an 
error with two additional ancilla qubits: a CNOT gate will duplicate the qubit in Z basis using the 
first ancilla qubit, the second ancilla will duplicate the qubit in X basis using two Hadamard gates 
and a CNDT gate (called an HCNOT gate), see Figure [2l[a). In this example, qubit 3 is the one affected 
by error in the last error correction cycle. 

2. Extended syndrome measurement. Measure the extended syndrome Sn+2 for an {n + 2) 

extended codeword consisting of the original codeword and the two extra ancilla qubit. Call E 
the result of the measurement of the extended syndrome Sn+2- When implementing this extended 
syndrome measurement, we need consider two aspects: (i) the entanglement between the additional 
ancilla and codewords should not be disturbed, since we need to disentangle the extra ancilla from 
the code and keep the codeword intact; and (ii) the new stabilizers (or generators) should also 
stabilize the codeword and satisfy all the stabilizers' requirements. In our approach, the additional 
ancilla qubit A is a copy of the control qubit in Z basis. To keep the entanglement, we copy all X 
operations on the control qubit to qubit A. Similarly, the additional ancilla qubit B is a copy of 
the control qubit in X basis. We replicate all Z operations performed on the control qubit to this 
qubit. Since it is in X basis, the Z operations will change into X operations. Consider the example 
of the Steane code and the time-correlated error on qubit 3, Figure [2l^b). In this case we replicate all 
X operations performed on qubit 3 to ancilla qubit A, and replicate all Z operations to additional 
ancilla qubit B as X operations. These operations will keep the entanglement intact during the 
syndrome measurement process. Also, it is easy to check that the new stabilizers (i) commute with 
each other, and (ii) stabilize the codewords. 

3. Disentangle and measure the two additional ancilla qubits. After the syndrome mea- 
surement, we use a circuit similar to the one in Figure [2l^a) to disentangle the two additional ancilla 
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Figure 2: (a) Duplicate the qubit affected error during the last cycle in both X and Z basis; (b) 
Extended syndrome measurement on both the codeword and the additional ancilla qubits. 

qubits from the codeword. The gates in this circuit are in reverse order as in Figure [^l^a). Since 
the entanglement is not affected during the extended syndrome measurement, this process will dis- 
entangle the extra ancilla from the codewords and keep the codeword in its original state. If qubit 
3 is affected by the time-correlated error, bit-flip, phase-flip or both, the disentanglement circuit 
will propagate these errors to the extra ancilla qubits: a CNOT gate will propagate a bit-flip of the 
control qubit to the target, and an HCNOT gate will propagate a phase-flip of the control qubit to 
the target qubit as a bit-flip. Figure [3l Therefore, measuring the additional ancilla qubits after the 
disentanglement will give us the information regarding the type of error qubit 3 has relapsed to. 



(a) 



4) — [K 



H • H 



H • H ^ Z 



(b) 



Figure 3: (a) A CNOT gate propagates a bit-flip of the control qubit to the target qubit (b) An HCNOT 
gate propagates a phase-flip on the control qubit to the target qubit as a bit-flip. 

4. Correct the errors according to the output of the extended syndrome measure- 
ment and of the measurement of the two additional ancilla qubits. Consider the following 
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scenarios: 



1. A "new" error affects a qubit i of the codeword, while qubit j,j ^ i was the one in error 
during the previous cycle. Then the measurement of the extended syndrome produces the 
same outcome as the measurement of the original syndrome and the new error is identified. 

2. The time-correlated error of qubit j occurs and there is no "new" error. The error propagates 
to the two ancilla qubits. The measurement of the original syndrome allows us to identify the 
single error. 

3. There are two errors, the time-correlated error of qubit j and a "new" error affecting qubits 
i,i ^ j of the codeword. Then the extended syndrome measurement gives the combination of 
these two errors: S = Ej © S^. The measurement of the two additional ancilla qubits reveals 
the type of the time-correlated error: 10 - bit-flip, 01 - phase-flip, 11 - bit-and-phase-flip, 
see Figure |31 Knowing that qubit j is in error and the corresponding error type we use the 
syndrome table to find the new error on qubit i. 

4. There is a single "new" error affecting one of the two extra ancilla qubits. This error may 
propagate to the qubit of the codeword the ancilla is entangled with. We treat this scenario as 
(1). For example, in Figure [21 a phase-flip of the ancilla qubit A propagates to qubit 3. The 
measurement of the original syndrome allows us to identify this error on qubit 3. 

5. There are two errors, the time-correlated error on qubit j of the codeword and an error on the 
additional ancilla entangled with qubit j. As we mention in scenario (3), the error on extra 
ancilla may back propagate to qubit j. The net result is that in the codeword, only the qubit 
j is affected by error: bit-flip, phase-flip or bit-and-phase-flip. This error is indicated by the 
syndrome and can be corrected as in scenario (1). 

The error correction process can be summarized as follows: 

1. If the syndrome S indicates an error on qubit j, the one affected by error during the last cycle 
(qubit 3 in our example), correct the error; 

2. Else, determine the outcome of a measurement of the two additional ancilla qubits: 

(a) If the outcome is 00, and E corresponds to a single error syndrome, correct this error; 

(b) If the outcome is either 01, or 10, or 11, qubit j in error during the previous cycle has 
relapsed. Then the syndrome is E^ = S (8 Ej. Knowing that qubit j is in error and the 
corresponding error type (10 - bit-flip, 01 - phase-flip, 11 - bit-and-phase-flip) we use 
the syndrome table to find the new error affecting qubit i. 

(c) Otherwise, there are two or more new errors in the system and they cannot be corrected. 



Example. Consider again the Steane sever-qubit code. Table [3] shows the syndromes for all single 
errors for the Steane 7-qubit quantum code. Assume that during the last error correction cycle qubit 
j = 3 was affected by an error. We use a CNOT gate to entangle qubit 3 with the additional ancilla 
qubit A in Z basis, and an HCNOT gate to entangle it with qubit B in X basis. Figure [2] 

1. E = 011000, indicates a phase-flip on qubit 3. We correct this error. In this case the measure- 
ment of the additional ancilla qubit could be non-zero, but this will not affect the codeword. 

2. E = 110000 and the outcome of the measurement of the two additional ancilla qubits is 00. In 
this case there is only one "new" error in the system, a phase-flip on qubit i = 6 as indicated 
by E. We correct this error. 

3. E = 110000 and the outcome of the measurement of the additional ancilla qubits is 01. In 
this case Ej = '^phase-fUp — = 011000. Therefore, the new error has syndrome E; = 
110000 ® 011000 = 101000; this indicates a phase-flip of qubit « = 5. Correct the phase-flips 
on qubit 3 and qubit 5. 
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Table 3: The syndromes for each of the three types of errors of each qubit of a codeword for the 
Steane 7-qubit code: Xi_7 bit-flip, Zi_7 phase-flip, and Yi_7 bit-and-phase flip. 



Error 


Syndrome 


Error 


Syndrome 


Error 


Syndrome 


Xi 


000001 


Zi 


001000 


Yi 


001001 


X2 


000010 


Z2 


010000 


Y2 


010010 


X3 


000011 


Z3 


011000 


Y3 


011011 


X4 


000100 


Zi 


100000 


Yi 


100100 


X5 


000101 


Z5 


101000 




101101 


Xe 


000110 


Ze 


110000 


^6 


110110 


X7 


000111 


Z7 


111000 


Yj 


mill 



The circuits described in our paper have been tested using the Stabihzer Circuit Simulator [5]; 
sample codes are provided in the Appendix. One can test different error scenarios using a code 
similar to the one in the Appendix. 

5 Summary 

Errors affecting the physical implementation of quantum circuits are not independent, often they 
are time-correlated. Based on the properties of time-correlated noise, we present an algorithm that 
allows the correction of one time-correlated error in addition to a new error. The algorithm can be 
applied to any stabilizer code when the two logical qubits | Ol) and | 1^} are entangled states of 2" 
basis states in 7i2'» . 

The algorithms can be applied to perfect as well as non-perfect quantum codes. The algorithm 
requires two additional ancilla qubits entangled with the qubit affected by an error during the 
previous error correction cycle. The alternative is to use a quantum error correcting code capable 
of correcting two errors in one error correction cycles; this alternative requires a considerably more 
complex quantum circuit for error correction. 
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APPENDIX 



# Sample code for testing the algorithm 

# Requires graphsim from 

# \protect\vrule wicithOpt\protect\href-Chttp: //homepage 

import graphsim 

import random random. seedO 

gr=graphsim . GraphRegister (22 , random. randint (0 , 1E6) ) 

# Prepare Steane I 0> codeword on qubit 0-6 
gr . hadamard (4) 

gr . hadamard C 5 ) 
gr. hadamard (6) 
gr.cnot (6, 3) 
gr . cnot (6 , 1) 
gr . cnot (6 , 0) 
gr . cnot (5 , 3) 
gr . cnot (5 , 2) 
gr . cnot (5 , 0) 
gr.cnot (4, 3) 
gr.cnot (4, 2) 
gr.cnot (4, 1) 

# Entangle qubit 3 with the extra ancilla 
gr.cnot (2,7) 

gr . hadamard (2) 
gr.cnot (2,8) 
gr . hadamard ( 2 ) 

# Errors happen on qubit ? 
gr.bitflip(2) 

gr .phasef lip(3) 

##### Begin current error correction cycle ##### 

# "Extended" Syndrome extraction 
for i in range (9, 15); 

gr . hadamard ( i ) 

# Z syndrome 
gr.cnot (9, 3) 
gr.cnot (9, 4) 
gr . cnot (9 , 5) 
gr.cnot (9,6) 

gr.cnotClO,!) 
gr.cnot (10, 2) 
gr.cnot (10, 5) 
gr.cnot (10, 6) 
gr.cnot(10,7) 

gr.cnot(ll,0) 
gr.cnot (11, 2) 
gr.cnot(ll,4) 
gr.cnotCll,6) 
gr .cnot(ll,7) 

# X syndrome 

gr . cphase (12 ,3) 
gr . cphase (12 ,4) 
gr . cphase (12 ,5) 
gr . cphase (12 ,6) 

gr . cphase (13, 1) 
gr. cphase (13, 2) 
gr. cphase (13, 5) 
gr. cphase (13, 6) 
gr.cnot(13,8) 

gr . cphase (14 , 0) 
gr . cphase (14 ,2) 
gr . cphase (14 ,4) 
gr . cphase (14,6) 
gr.cnot (14, 8) 



# "Disentangle" the extra ancilla 
gr . hadamard ( 2 ) 

uibk . ac . at\gti:iai5t(fEQBjl3/work/graphsim .htmlXhttp : //homepage . uibk . ac . at$\sim$c705213/work/g 
gr . hadamard ( 2 ) 
gr.cnot (2, 7) 

# Measure and output the extra ancilla 

print gr .measure (7) 
print gr .measure (8) 

#gr .print_adj_list () 
#gr . pr int_stabilizer ( ) 
print 

# Measure and output syndrome 
for i in range (9, 15): 

gr. hadamard (i) 
print gr.measure(i) 
print 

############# error correction ############### 

# Error correction according to the syndrome and extra ancilla 
gr .bitflip(2) 

gr .phasef lip (3) 

##### Normal syndrome extraction to double-check ##### 

for i in range (9+6 , 15+6) : 
gr .hadamard(i) 

# Z syndrome 
gr.cnot (9+6, 3) 
gr.cnot (9+6,4) 
gr.cnot (9+6,5) 
gr.cnot (9+6, 6) 

gr.cnot (10+6,1) 
gr.cnot (10+6, 2) 
gr.cnot (10+6, 5) 

gr.cnot(10+6,6) 

gr.cnot(ll+6,0) 
gr.cnot(ll+6,2) 
gr.cnot(ll+6,4) 
gr.cnot(ll+6,6) 

# X syndrome 

gr .cphase(12+6,3) 
gr.cphase(12+6,4) 
gr . cphase (12+6 ,5) 
gr.cphase(12+6,6) 

gr.cphase(13+6,l) 
gr.cphase(13+6,2) 
gr.cphase(13+6,5) 
gr.cphase(13+6,6) 

gr.cphase(14+6,0) 
gr.cphase(14+6,2) 
gr . cphas e ( 14+6 , 4 ) 
gr. cphase (14+6, 6) 

# Measure and output the syndrome 
for i in range (9+6, 15+6) : 

gr .hadamard(i) 
print gr .measure(i) 
print 
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